Separation of energy scales in the kagome antiferromagnet TmAgGe: a 
magnetic-field-orientation study up to 55 T 
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TmAgGe is an antiferromagnet in which the spins are confined to distorted kagome-like planes 
at low temperatures. We report angle-dependent measurements of the magnetization M in fields of 
up to 55 T that show that there are two distinct and separate energy scales present in TmAgGe, 
each responsible for a set of step-like metamagnetic transitions; weak exchange interactions and 
strong crystalline electric field (CEF) interactions. Simulations of M using a three-dimensional, 
free-energy minimization technique allow us to specify for the first time the physical origin of the 
metamagnetic transitions in low, in-plane fields. We also show that the transitions observed with 
the field perpendicular to the kagome planes are associated with the CEF-split multiplet of Tm. 

PACS numbers: 75.25.-(-x, 75.10.-b, 75.30.Gw, 75.30.Kz 



Interest in antiferromagnetic (AF) systems based on 
structures with possible frustration, such as triangu- 
lar, kagome or pyrocholore lattices, has recently bur- 
geoned, as such materials often exhibit novel cooperative 
phases One example, exhibiting complex metamag- 
netic behavior, is TmAgGe, a member of the i?AgGe 
(i?=Tb, Dy, Ho, Er, Tm and Y) family of metallic com- 
pounds. TmAgGe adopts a layered, distorted-kaeome 
lattice similar to the ZnNiAl structure (Fig. [T|) [3, 3- 
In this paper, we show that a model based on a six- 
spin repeating structure can account for the majority of 
magnetic data on TmAgGe, including the field {H) posi- 
tions of the metamagnetic transitions and the magnitude 
of the magnetization Af , when H is applied within the 
kagome planes. Though TmAgGe exhibits few of the 
features conventionally associated with frustration, the 
model shows that it conforms to the fundamental defini- 
tion of a frustrated system [l| : the geometry of the lattice 
precludes the simultaneous minimization of all of the in- 
teraction energies, in this case antiferromagnetic (AF) 
next-nearest-neighbor (NNN) and ferromagnetic (FM) 
nearest neighbour (NN) exchange with a characteristic 
energy scale ~ 4 K. We also report the first observation 
of a series of high-field metamagnetic transitions when 
H is approximately parallel to c; these are due to field- 
induced level crossing within the crystalline-electric-field 
(CEF) split Tm3+ J = 6 (4fi2, SHg) multiplet. Here, the 
energy scale is ~ 100 K. Based on these two energy scales 
(i.e. exchange ~ 4 K, CEF ~ 100 K), we can produce a 
complete quantitative phase diagram for TmAgGe. 

Much of the essential physics is captured in Fig. [2l 
Fig- 121 a-) shows the metamagnetic transitions (peaks in Xi 
corresponding to steep rises in M) at several values of 4>, 
the angle between H and the c-axis. These low- if tran- 
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FIG. 1: Projection of the distorted kagome structure of 
TmAgGe in the a&-plane. Tm: red, Ag: silver, Ge: gold. The 
model in this paper has spins 1-6 as its repeating structure. 
Nearest neighbors 1, 2 & 3 and 4, 5 & 6 have AF interactions; 
3 has FM interactions with 4 and 5; because every second tri- 
angle of spins is considered equivalent, 2 has FM interactions 
with 4 and 6. No FM interactions act on 1 in this model. 

sitions, determined by the ~ 4 K energy scale, depend 
only on the component of H parallel to the a5-planes. 
At higher H , the new series of transitions is observed 
(Fig. mb)); as shown below in the discussion of Fig. |4l 
these result in M approaching its full saturated value of 
per Tm'^"'" ion. These transitions, determined by the 
^ 100 K energy scale, depend only on the component 
of H perpendicular to the a6-planes. The (iJ, 0) phase 
diagram is in Fig. [2Ic); note that the low- and high-iJ 
transitions can be seen in the same field sweep, showing 
that they are distinct phenomena, separated at all (j). 

Oriented single crystals of TmAgGe (for growth de- 
tails, see Q), ~ 0.7 x 0.7 x 2 mm^, were used in 
compensated-coil M measurements in a 65 T pulsed mag- 
net at NHMFL. In one probe, the sample can be inserted 
into and extracted from the coil in situ, enabling a mea- 
surement of the sample's M . In the second probe, the 
sample is in a coil that tilts so that the angle-dependence 
of X = AM/ AH can be recorded. This was used to ob- 
tain the data in Fig. [51 The inclination is deduced by 
comparing the voltage induced in an empty coil on the 
tilting platform, and that from another empty, static coil. 
Both probes were placed in ■^He cryostats (T > 500 mK). 
Pulsed-field data were calibrated against results from 
Quantum Design MPMS systems {ixoH < 7 T; T > 2 K). 

The sharp features in x observed in Fig. [Ifa) corre- 
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spond to steep increases of M on either side of rounded 
plateaux. Examples are shown in Fig. [2Id) for (f) — 90° 
and 6* = (i?||[110]) and 9 = 30° (i/||[120]); here, 6 is 
an azimuthal angle coordinate, for H rotating in the ab 
plane. In the former data, a broadened plateau is centred 
on M « 2.3/iBTm^^ and in the latter, a smaller, rounded 
shelf can be discerned at M « 2.0/XBTm^^. These data, 
and others recorded at fixed 9, reproduce quantitatively 
the M{H,9) measurements in Fig. 24 of Ref. 3. To 
model such M data, and to predict the {H, 9, (j)) phase 
diagram of TmAgGe, we use the following Hamiltonian: 
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•k — Jfi^p ^ Si ■ Sj + JpM ^ Si • Sj 

+A^(5f)2-g^BES.-B. (1) 

i i 

Here, Jaf > and Jfm < are NNN AF and NN 
FM exchange constants, is the component of the 
ith spin along the local easy axis 0] and A < is 
the anisotropy energy. The 6-spin repeating structure 
(Fig. [T|) is sufficient to reproduce the majority of the fea- 
tures of M(i7, 61). 

Having established the interactions and the number of 
spins involved, the simulation involves using the down- 
hill simplex method to find the minimum free-energy 
spin configuration for a particular value and direction of 
H. The parameters Jaf, Jfm and A are then adjusted 
to quantitatively match AI{H,9,(p) data {e.g. Figs.HJd) 
andini^b)), a process that provides a tight constraint of 
the values. Note that in contrast to previous phenomcno- 
logical geometrical models of TmAgGe Q, there is no 
need to impose a starting spin configuration; our model 
automatically finds a reasonable arrangement. 

Simulations are shown with data in Fig.[2Id), and cor- 
responding spin configurations are given in the (iJ, 6) 
phase diagram {(j> = 90°) of Fig. [2Ia) for the various 
states. In the AF groundstatc, the spins on each trian- 
gular plaqucttc lie in the non-colinear directions dictated 
by the anisotropy, with spins pointing inwards and out- 
wards alternately on adjacent triangles in order to mini- 
mize the energy of the FM term in the Hamiltonian. As 
H increases along [110] {9 — 0) the spins remain parallel 
to their easy axes, but three of the spins flip towards the 
direction of H, while two flip back to minimize the AF 
term's energy. This forms the M2 state Q , corresponding 
to the plateau in M in the model and the broadened step 
in the data. As H increases further, the AF interactions 
are overcome and the spins rotate towards H, leading 
to the crystal-field-limited saturated paramagnetic (CL- 
SPM) state. For i/||[120] {9 = 30°), the M2 state again 
forms (plateau in model, inflexion in data) but the CL- 
SPM state is not realized, as one of the spins is perpen- 
dicular to, and so cannot couple with, H. Instead, the 
FM interactions force the moments into the M3 state. 

In addition to the absolute size of M within the various 
states (color scale). Fig. |31[a) shows the 6'-dependences of 
the critical fields marking the boundaries between the 
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FIG. 2: (a),(b) Magnetic suceptibility x = dM/dH of 
TmAgGe (T « 500 mK) measured in pulsed magnetic fields 
at various values of the inclination angle (j). Peaks in x cor- 
respond to rapid rises of M (transitions) between broadened 
steps. Data are offset for clarity, (c) Critical fields of the 
magnetic transitions (T « 500 mK) plotted as a function of 
<f). AF — antiferromagnet; CL-SPM — crystal-field limited 
saturated paramagnet; SPM — saturated paramagnet. Dot- 
ted fines are best fits to l/|sin0| (low fields) and l/lcos^ 
(high fields), (d) M predicted by the model (see text) with 
H along the [110] and [120] directions compared with experi- 
mental pulsed-field data. 

spin configurations (thick lines). These are obtained by 
equating the free energies of the states, yielding: 



oAF^Ma 
r>M2-»M3 
oMa^CL-SPM 
^M3-»CL-SPM 



/3 Jaf/ cos 6*; 

/3(2Jaf- Jfm)/(2cos(0- V3)) 
/3(Jaf - Jfm)/ cos 0; 
-/3JFM/(2cos(6' + 7r/3)). 



(2) 



Here, /3 — J/gj^B, where J = 6 and gj — 7/6 are re- 
spectively the total angular momentum quantum number 
and Lande g- factor for Trn'^^. 

Eqs. [5] are similar to those extracted from M{H, 9) 
data in Ref. 3], but with explicit prefactors; comparison 
of Eqs. with the experimental critical flelds observed 
in the current study and in Ref. 3] forms a tight con- 
straint on the model parameters Jaf and Jfm- We find 
Jaf = 0.054±0.003 K and Jfm = -0.064±0.003 K. The 
energy ranges of the first two terms in Eq. ([T]) are ± J^ Jaf 
and ± J^ Jfm; thus the energy scale at which the exchange 
interactions are important is 4 K. Not only does this 
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energy scale account for the sizes of the critical fields of 
the low-field metamagnetic transitions (Figs.[2I^a),(d)); it 
is also responsible for the H = transition at Tm ~ 4.2 K 
into the state that we have labelled. Although the AF 
state has zero moment, the model shows that the FM 
interactions are ultimately responsible for its spin config- 
uration. This, and the plethora of other states observed 
in the data (Figs. EJa), (b)) result from competition of 
FM and AF interactions; the geometry of the lattice pre- 
cludes the simultaneous minimization of all of the inter- 
action energies. In this respect, TmAgGe conforms to the 
most fundamental definition of a frustrated system [l| . 

Whilst our 6-spin system (Fig. [1]) and Hamiltonian 
(Eq. [T|) describe most of the low-field M{H, 9, (p) data 
presented in this paper and in Ref. Q , it is necessary to 
make two comments. First, the rise in M observed be- 
tween the low-field AFM state and the broadened step 
associated with the M2 spin configuration (Fig. ISl^d)) is 
attributed elsewhere [Sj to a state labelled Mi; our model 
does not predict Mi, instead showing a direct AFM-M2 
transition (fig. [Sja)). Second, to achieve a quantitative 
reproduction of all the other states observed, it is neces- 
sary to suppress FM interactions between spin 1 in Fig.[T] 
and its NNs: if these interactions are switched on in the 
simulations then the angular region over which the M3 
state is observed is severely reduced and the M2 state 
disappears completely, in contradiction to experiments. 
The most likely explanation is that our 6-spin model is 
a subset of the actual repeating magnetic structure in 
TmAgGe. Support for this view comes from the struc- 
tures inferred from analysis of M{H, 9) data in Ref. 
which suggest that > 18 spins are required to reproduce 
M within the Mi state. Nevertheless, we emphasise that 
the 6-spin repeating system (Fig [U Eq. [T]) is able to de- 
scribe almost all of the data, especially those at high 
H, and represents a tractable model for exploring spin 
physics in layered kagome systems of this kind. 

We now turn to the other energy scale in TmAgGe, 
associated with the CEF interactions. Though the mo- 
ments are strongly inclined to lie along the easy axes, the 
finite A means that they will cant towards H, leading to 
the gradual increase in M as H rises seen in the CL- 
SPM state (Fig. [2l^d)). Using this gradient, we find that 
A = —4.6 ±0.1 K. This explains the large negative Curie 
temperature for H\\c (0c); the Zeeman and anisotropy 
terms dominate in Eq. [T]for H\\c, and in this limit one 
can show that the high-T 1/x is linear with a T-axis 
intercept of 9^ « -A(2J + 3)(2J - 1)/15 « -50 K in 
reasonable agreement with experiment (0c ~ —76 K [4']). 

The third term of Eq. [1] has a magnitude |^^A| ^ 
170 K, forcing a non-colinear or "compromise" structure 
on the Tm moments at Ts far in excess of those at which 
the AF interactions become important. This compromise 
structure, in which the moments lie in the ab planes 120° 
apart, could potentially possess a number of different de- 
generate spin configurations [l|. However, in TmAgGe 
the model shows that the degeneracy is lifted by the FM 
NN interactions, leading to the states at low Ts and the 
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FIG. 3: {H, e) (a) and [H, if)) (b) phase diagrams of TmAgGe 
predicted by the model (see text). The color keys show \M\. 
(a) Spin configurations are shown for each state; solid lines 
show critical-field 6'-dependences from Eq. [2l The diagram 
agrees well with data in this work (e.g. Fig. [2jd)) and in 
Fig. 25 of [3. Note that the diagram has a 9 periodicity of 
60°. (b) Points are experimental critical fields (Figs. [2ja), 
(c)); dotted lines are fits of data to 1/| sm(j>\ {9 — 20°). 

magnitude and sign of Qab ~ -t-7.5 K Q. As mentioned 
above, these states are then accessible using small in- 
plane {4> = 90°) Hs, accounting for the metamagnetic 
transitions (Figs. [IJd) and [3]) |71]. 

Fig. [3fb) shows the {H, 0) phase diagram of TmAgGe 
simulated using the above parameters, together with ex- 
perimental data. In agreement with experiment, the rela- 
tively large [Aj confines the spins to the ab planes, so that 
the \ow-H, metamagnetic transitions depend only on the 
component of H in the planes; the critical fields scale as 
1/1 sin 01 0. Note that the simulation uses 9 = 20°, in 
agreement with the nominal orientation of the sample. 

We now turn to the high-i/ metamagnetic transi- 
tions observed when H is applied out of the ab planes 
(Fig. [2Ib)). As mentioned before, peaks in x corre- 
spond to steep rises in M; this is shown in Fig. HJ^a) 
for H\\c{9 = 0). Eq. (P) predicts that for H \\ c, the 
low-T M should increase linearly in H with a gradient 
given by (77/13/2 A, reaching its saturated value of 7fiB 
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FIG. 4: (a) M{H,T) of TmAgGe {H \\ c) showing the 
transitions attributed to CEF-spht energy-level crossing. In- 
set: critical fields versus T. (b) Calculated M with H \\ c 
found by diagonalizing H\ (see text), and experimental data 
(T — 500 mK for both). The Stevens parameters were B2 = 
1300 mK; Bl = -3.1 mK; B| = 19 mK; and = 0.0068 mK, 
all the rest being zero. Inset: calculated 1/x using the same 
at two different directions of H, compared with data ob- 
tained at 1 T 

per Tm'^+ ion at ^qH « 70 T. The experimental M in 
Fig. [3)Ja) does increase linearly at low H (albeit at low 
T with a smaller gradient than predicted); however, at 
/iQ-ff w 35 T the steep rise in M mentioned above is ob- 
served, and M approaches its saturated value. This is 
similar to what has been observed in high-i/ measure- 
ments of HoNizBzC For T < 20 K, there is a two- 
or three-fold structure within the rise in M; at higher T, 
these structures are not resolved, but a single transition 
in M is visible to T « 70 K (Fig. IDJa), inset). Ts of 20 K 
and 70 K are several times higher than Tm ~ 4.2 K, indi- 
cating that the transition(s) in M in Fig. |4]cannot easily 
be attributed to exchange interactions. Indeed, as one 
tilts H away from c (Figs. [DJa), (b)), the low-i? mctani- 
agnetic transitions that can be confidently attributed to 
the competing FM and AF exchange interactions are ob- 
served within the same field sweeps as the transition(s) at 



^lqH « 35 T; this strongly suggests that the mechanisms 
for the low- and high-i? transition(s) are distinct. 

The field of /io-ff w 35 T suggests an energy scale for 
the high-_ff transitions similar to the anisotropy term 
in Eq. [T] The likely candidate is CEF splitting of the 
J = 6 multiplct, which causes the easy-axis anisotropy 
seen at low T. As stated above, such splittings would 
be of the order of | J^A] - 170 K (Eq. [H third term). 
When the Zeeman energy = gj^^JB becomes of this 
order, level crossings are expected; this occurs when 
fioH « 36 T, in agreement with the data. The Tm 
ions are located at sites with orthorhombic symmetry 
C2v (2mm) i and so we consider the Hamiltonian 
Hi = E„=2,4,6 E:=o - 5Mb S. • B, a sum of 

CEF [llj and Zeeman terms, where the and are 
the CEF parameters and Stevens operators, respectively. 
The O™ with m 7^ contain the angular momentum 
raising and lowering operators; thus, if the correspond- 
ing CEF parameters are large enough, these terms will 
lead to substantial mixing of the Jz energy levels. 

By diagonalizing Hi and using expressions for x and 
M from Ref. [l^, we find values for the Stevens param- 
eters such that the low-T M with H \\ c is a reasonable 
match to the experimental data (Fig. IHb)). The low-_ff 
c-axis M is well reproduced, as is the transition to satu- 
ration around 35 T. However the three steps within this 
transition are not found in the calculations. The inset 
in Fig. (Ub) shows the calculated low-field 1/x versus T. 
Whilst the calculations reproduce the inset data quali- 
tatively, the quantitative agreement is less good than in 
Fig. mb), partly because at low T the calculated 1/x is 
very sensitive to (fi; the same is true to a lesser extent of 
the experimental data. Note that the values in the 
caption to Fig. [4] merely represent a possible solution to 
the problem; our intent here is to show that CEF split- 
ting of the J = 6 multiplet alone can account for the 
primary features of the T dependence of the c-axis x f^nd 
M {H). A further constraint on the i?'" parameters must 
await inelastic neutron scattering experiments. 
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